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PREFACE 


In order to accurately predict the state (i.e. , position and 
velocity) of a near-Earth satellites the equations governing the 
satellite's motion must be integrated numerically because the pertur- 
bations produced by atmospheric drag, an oblate Earth and celestial 
bodies preclude the finding of a closed form solution. Generally, the 
integrated states differ from the actual states, and these differences 
(state deviations or enrors) increase with time. These state devia- 
tions are influenced to a high degree by the coordinate system in 
which the equations of motion are expressed. Since the tiro-body 
influence dominates the satellite's motion and since a closed form 
snii.rHf'n to the two-bodv DToblem is available, an analysis of how two- 
body motion is “affected by different coordinate systems would provide 
some indication as to how the more complicated perturbed motion is 
affected by coordinate system choice. 

The purpose of this investigation is to determine the 
influence which various coordinate systems have upon the error in 
the numerically determined value of the state for tiro-body motion 
using several integration step sizes and orbit shapes. This is accom- 
plished by numerically integrating the equations of motion for each 
coordinate system under investigation. At specified time intervals 
an error check is made by transforming the respective position and 
velocity vectors, to a reference coordinate system. Then a comparison, 
of the integrated state and a reference state generated using the 

iii 



solution of the two-body equation is made. In this manner the 
relative state errors for each system can be obtained, and the res- 
pective error norms computed. After varying the integration step 
size and the orbit shape, the error norms are analyzed m order to 
determine the sensitivity of the relative error magnitudes to step 
size, orbit shape and coordinate system choice. 

The author wishes to express his sincere appreciation to 
Dr. W. T. Fowler for bringing this area of study to 1 the author's 
attention, for his help and guidance throughout the course of this 
study and for serving as the supervising professor. The author 
wishes to thank Dr. B. D. Tapley and Dr. D. G. Hull for serving on 
the supervising committee, Dr. G. J. Lastman for his valuable sug- 
gestions and advice, Mr. Jack Hart for initial programming assistance 
and Mr. Jim Wiley for checking the equations of motion. Finally, the 
author wishes to thank his parents for their help and encouragement 

throughout his education and particularly his wife, Janice, for her 
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love, patience and understanding. 

The author is pleased to acknowledge the support of the 
National Aeronautics and Space Administration under Grant NSG-5S1. 
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CHAPTER I 


INTRODUCTION 


1.1 Purpose of Investigation 

Since the launching of the first artificial Earth satellite 
in 1957, man has been concerned with the problems encountered in the 
computation of such satellite orbits. In general, this is a very 
difficult task to perform accurately over long time intervals. This 
task is usually carried out numerically because closed form solutions 
do not exist to the complete ordinary second-order nonlinear differen- 
tial equations which govern the satellite's motion. In order to 
describe this motion as accurately as possible, the equations of 
motion must include terms accounting for all forces influencing the 

f , . *..* . T T - „ , * ► * - * V f ~ 

C*- W-L W AVh> Ol °* L,A 5 

vehicles are non thrusting, their mass is much less than that of the 
Earth, and their orbits are such that the effects of the perturbing 
force effects are small, a fairly accurate orbit description can be 
obtained from two-body considerations. By assuming an inverse 
square force field, the two-body equation of motion can be expressed 
as 


d r -yr 

— 7 — r 

dr r 


CD 


where y is the Earth's gravitational constant, r is the position 
vector measured from the center of the Eartli and r is the magnitude 
of the position vector. A closed form so]ution to Equation (1) is 
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2 


given by 

r = a(l-e 2 ) 

1 + e cos <{> 


( 2 ) 


where a is the semimajor axis, e is the eccentricity and d> is 
the true anomaly [see Figure 1) . Equation (2J is the equation of a 
conic section in polar form and consequently represents a circle 
(e - U) , an ellipse (0 < e < 1) , a parabola (e =1) or a hyperbola 
(e > 1] . Since closed orbits [e < 1) ate of the most interest, only 
the orbit geometry associated with e < 1 is shown in Figure 1. Also 
shown in Figure 1 is the eccentric anomaly, E, which related the 
elliptic motion to the equivalent circular motion. Hie eccentric 
anomaly is related to the true anomaly, <j>, by 


tan E 


(1-e 2 )" 2 sin 
1 + e cos <j> 


(3) 


The eccentric anomaly is related to the time in orbit, t, from 
perigee passage, t^, by 


t - tp = (— ) (E-e sin E) 


Whenever the states for a two -body orbit are obtained by 
--integrating the equations of motion numerically, the states obtained 
differ from those obtained via the two -body solution. Hie error 
which accumulate in the integrated states. is a function of the type 
of orbit being integrated, the coordinate system in which the equa- 
tions of motion are formulated and the size of the integration time 




Figure 1. Elliptic Orbit Geometry 
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step. The manner in which the above three parameters are combined 
affects both the accuracy and the computational time (which can be 
exceedingly long in some cases such as orbital^lifetime and epher- 
meris predictions) , Because of the expense involved in computer 
operation, considerable emphasis is placed on efficient programming 
in order to reduce running time, thus conserving cost, while still 
achieving the greatest accuracy possible. 

The purpose of this study is to investigate the relative 
error propagation in the state variables resulting from the numeri- 
cal integration of near-Earth two-body orbits. The influence of 
coordinate system, integration step size, and orbit shape on the 
computational efficiency is examined with the emphasis being placed 
on accuracy and reduced running times. 

I . 2 Previous Studies 

1 * 

A previous study has been conducted by Gerber and Lewallen 
in which the generation of relative state error magnitudes as a 
function of orbit shape and integration step size for rectangular 
and spherical coordinate systems was investigated. ' 

This study first determined the effects of integration step 
size on circular orbits. The step sizes considered ranged from 15 
seconds to 120 seconds. From these results a step size of 30 seconds 
was selected for use in studying eccentricity variations which ranged 

* Numbers appearing m the text as superscripts indicate 
references listed in the Bibliography. 
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for 0 to 0.3. The orbits were integrated over an interval of only 
60,000 seconds because of computer core storage limitations. The 
computer used in this study was a Univac 1108 with a 36 -bit word 
length. Tills study was carried out in single-precision arithmetic 
with partial double precision in the integrator only. Other 


assumptions and simplifications were essentially the same as the 
ones which will be used in the present investigation and will be 
discussed later. 

The preliminary conclusions reached by Gerber and Lewallen 
were as follows. First, the velocity error norms and the position 
error norms exhibited the same types of behavior for both step size 
and eccentricity variations. Second, the errors in the spherical 




:11 beh 


1 *t noo-rlir t.ri 'f'Vt 


accumulated integration time while the errors in -the rectangular 
coordinate system exhibited a fairly constant amplitude periodic 
variation superimposed on an exponential type increase. The rectan- 
gular coordinate errors were generally larger than the spherical 
coordinate errors and the differences increased as the orbital 
eccentricity increased. Typical terminal time results are given in 
Table 1 where the quantity A/P refers to the ratio of apogee height 
to perigee height with height being defined as the distance above 
the Earth's surface and the quantity ISS refers to the integration 
step size. Third, the results indicated that a "best" numerical 
integration step size (a compromise between round-off and truncation 
errors) could exist for each coordinate system. The "best" step size 



TYPICAL RESULTS - 

GERBER AND LEIVALLEN 

Error in Position at 60,000 Seconds (feet) 

Coordinate System 

A/P = 1 

A/P = 2 

Rectangular 

1250 

960 

Spherical 

140 

265 


I Prnl^l pm PuromAforQ I 

) * ‘ ' l 


ISS = 30 seconds ft = 0.1 radian 

i = 0.5 radian oj = 0.1 radian 


TABLE 1 
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appeared to be larger for the spherical coordinates than for rectan- 
gular coordinates. Finally, it became apparent that the errors in 
the state were computer word length dependent, thereby suggesting 
that a comparison of results obtained using full double precision 
and/or other computers be made. 

2 

Later unpublished work by Rainbolt and Lewallen included 
analysis of the circular cylinder, parabolic cylinder, paraboloidal 
and elliptic cylinder coordinates as well as the rectangular and 
spherical coordinates. The elliptic cylinder coordinates used were 
not rectified to the orbit plane (i.e., did not take advantage of 
the particular form and properties of the coordinate system) . This 
work was conducted on the Univac 1108 computer using full double- 


0,1 ATI ovThhvnati <-• I)rp1 4-vAtn 

investigators indicated that the parabolic cylinder coordinates ex- 
hibited smaller error generation tendencies than any of the other 
coordinate systems studied. Typical terminal time results for this 
investigation are given in Table 2. 


1.3 Scope of Investigation 

In the present study, five coordinate systems (rectangular, 
spherical, circular cylinder, and two forms of elliptic cylinder 
coordinates) are compared in order to determine their efficiency for 
numerical integration of orbital equations. These five coordinate 
systems were chosen from numerous coordinate systems considered in a 
preliminary screening process. Various combinations of eccentricity 
and step size are investigated in order to determine the sensitivities 



TYPICAL RESULTS - 

RAINBOLT ANDlLEWALLEN 

Error in Position at 90,000 Seconds (feet) j 

Coordinate System 

A/P = 1 

A/P = 2 

Rectangular 

22.4 

11.0 

Spherical 

21.8 

20.5 

Paraboloidal 

17.0 

14.5 

Elliptic Cylinder 

17.5 

12.5 

Circular Cylinder 

13.0 

8.0 

Parabolic Cylinder 
. 

. . _ . 

5.8 



8,0 

1 

Problem Parameters 


ISS = 15 seconds 

& = 1.0 radian 

i = 0.5 radian 

co = 1.0 radian 


TABLE 2 




of the coordinate systems studied to orbit shape and integration 
step size. The eccentricities range from 0 to 0.2 while the in- 
tegration step sizes range from 15 seconds to 120 seconds. 



CHAPTER II 


ANALYSIS 


11. 1 Simplifying Assumptions 

Tiie assumptions made concerning the orbits treated in this 
study are as follows. The types of orbits considered are high 
enough so that the atmospheric drag can be neglected but low enough 
so that the celestial body perturbations can also be neglected. The 
orbiting satellite is nonthrusting and it is assumed that the Earth 
can be represented by a spherical body of homogeneous mass distribu- 
tion. 

11. 2 Selection of Variables and Constants 

In- order to assure that realistic orbit shapes and integra- 
tion step sizes are used in this investigation, care must be taken 
in selecting the values of these parameters. The orbits chosen 
should be representative of two-body motion, and the time stpes 
should be realistic ones. 

The motion of a body moving in a two-body orbit can be 
defined by specifying six quantities (orbital elements) . The stan- 
dard orbital elements are the semimajor axis, a, the orbital eccen- 
tricity, e, the orbit inclination, i, the longitude of the ascending 
node, R, the argument of perigee, w, and the time of perigee passage, 

t . 

P 

The function of the first two elements, a and e, is to 


10 




11 


define the orbit's shape. For this study a and e are obtained by 
specifying apogee and perigee for the orbit under consideration 
(see Table 3) . 

The next three elements, i, fi, and id locate the orbit with 
respect to the inertial reference coordinate system and remain con- 
stant throughout the course of this study. These variables have 
been arbitrarily assigned the following values: i = 0.5 radian, 

P = 0.1 radian and id = 0.1 radian. The remaining element, t , is 
set equal to zero since this element has no special significance 
for the present study as it is just a time reference. 

Remaining to be specified is the integration step size, 

ISS. The values considered for this parameter are 15 seconds, 30 
ScCoiiuS , 00 ScCuliOS dUu. iiiese v'H.aucs rtlfc tiMcd „i ui' 

all orbits studied. The relationship between integration step size 
and the associated error checks is given in Table 4. 

Other quantities which must be specified also are the total 
integration time, IT, the Earth's gravitational constant, p, and 
the Earth's radius, R . These values along, with the constant 
orbital elements are given in Table 5. 

II. 3 Selection of Coordinate Systems 

Prior to selecting a set of coordinate systems for detailed 
analysis, numerous coordinate systems were considered. Among those 
considered were rectangular coordinates, two foims of spherical 
coordinates, two forms of circular cylinder coordinates, (henceforth 



VARIABLES 

Perigee Ht. 
(miles) 

Apogee Ht. 
(miles) 

A/P 

Eccentricity 

Semimajor Axis 
(feet) 

100 

100 

1 

0.0 

2.145E+7 

100 

200 

2 

0.0125 

2.172E+7 

100 

500 

5 

0.0469 

2.2S1E+7 

100 

1000 

10 

0.0997 

2. 383E-* 7 

100 

2000 

20 

0.1895 

2.649E+7 


ISS = 15, 30, 60, 120 seconds 


TABLE 3 


STEP SIZE - ERROR CHECK RELATIONS 


Steps 


Step Size 
(seconds) 

Number of 
Integration Steps 

Number of 
Error Checks 

IS 

12000 

300 

30 

6000 

300 

60 

3000 

300 

120 

1500 

300 


Between Between Checks 



Integration Time = 180000 seconds 


TABLE 4 



















CONSTANTS 


i = 0.5 radian 
= 0.1 radian 
w = 0.1 radian 


t =0.0 seconds 
P 

y = 1.407653 x 10 16 £t 3 /sec 2 

R = 3963.115 miles 
e 
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called circular coordinates) , three forms of elliptic cylinder 
coordinates, (henceforth called elliptic coordinates), ellipsoidal 
coordinates, and conical coordinates. 

In the following paragraphs each coordinate system will 
be- briefly discussed, and the reasons for either rejecting or re- 
taining the system will be given. In an effort to increase com- 
putational efficiency, any substitution or coordinate system 
orientation which significantly simplifies either the coordinate 
transformations or the equations of motion is employed, 

I I. 5.1 Reference Coordinate System 

Hie coordinate system chosen as a reference system in 
this study is an inertial Earth- centered rectangular coordinate 
system with the XY plane lying in the equatorial plane. This 
coordinate system, shown in Figure 2 , is labeled the (X,Y,Z) sys- 
tem. 


11. 3. 2 Rectangular Coordinates 

This coordinate system, also designated by (X,Y,Z), is 
in common use and is therefore retained for the purposes of this 
study. The rectangular Cartesian system used coincides with the 
inertial reference coordinate system (See Figure 2) . 

11. 3. 3 Spherical Coordinates 

Two forms of spherical coordinates (Figure 3) were con- 
sidered. They were the standard spherical coordinates (p ,o , t) 
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and a nonstandard set of coordinates (p,cos o,cos r) . Although the 
equations of motion for the nonstandard system were simple, an 
ambiguity in the sign of the sine term exists. In order to resolve 
this sign ambiguity, a sign convention employing a series of checks 
on the coordinate values was developed. A preliminary computational 
check indicated that integration of the system containing the angles 
was inherently more accurate than integration of the system contain- 
ing the cosines of the angles. Due to the accuracy difference 
(probably due to the oscillatory nature of the cosine) the nonstandard 
system was discarded while the standard set was retained. Again, in 
order to eliminate a coordinate transformation, the spherical coor- 
dinate axes (X,Y,Z) of Figure 3 were aligned with (X,Y,Z) axes of the 


11.3.4 Circular Coordinates 

Two forms of the circular coordinates were considered 
(Figure 4) „ They were the standard set r, G , z and a nonstandard 
set (r,cos e,z). Hie nonstandard set exhibited accuracy problems 
equivalent to those of the nonstandard spherical coordinates and 
was rejected. The standard circular coordinates were retained and 
as in the previous cases, the coordinate axes (X,Y,Z) shown in 
Figure 4 were aligned with (X,Y,Z) axes of the reference system. 

11. 3. 5 Elliptic Coordinates 

Three sets of elliptic coordinates (Figure 5) were con- 
sidered. Set 1 designated by (ccpc^,^) and set 2 designated by 




(G;L>$2>&3^ werc rectified to the orbit plane while set 3 designated 
by ^ y i ,y 2 jY 3 -) was re f erence< i to the inertial reference system. Hie 
coordinates of set 1 are defined by 

= d cosh u (*2 = v aj = z 

where d is a constant related to the coordinate geometry. The 
coordinates of set 1 are measured in the orbital plane. Thus, the 
elliptic axes (X,Y,Z) shown in Figure 5 were aligned with the 
x,y,z axes of Figure 2. The coordinates of set 2, also measured in 
the orbit plane, are defined by 

~ d cosh u $2 “ cos v 63 = z 
The coordinates of set 3 are defined by 

Y 1 = u Y 2 = v Y 3 = z 

where (u,v,z) are the standard elliptic coordinates shown in Figure 
5 and are measured in the equatorial (XY) plane. Thus, the (X,Y,Z) 
axes for this set are aligned with the (X,Y,Z) axes of Figure 2. 

Note that in order to obtain the (X,Y,Z) (reference system) components 
of coordinate sets 1 and 2 , the (x,y,z) components in the orbital 
plane (see Figure 2) must first be obtained from the respective 
elliptic coordinates and then transformed to the reference coordinate 
system using the transformations given in Appendix A. 5. 
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The standard procedure for determining the initial values 
of the elliptic coordinates is to solve Equation (8) for the par- 
ticular values 

X = d cosh u cos u Y = d sinh u sin v Z = 3 (8) 


where u, v and z are the elliptic coordinates as they are ordinarily 
defined (set 3] and d is a constant related to the coordinate 
geometry (see Figure 5) . Although the above transformations are 
straightforward, the inverse transformations (necessary in set 3) 
are quite complicated. However, it will be shown in the following 
paragraphs that a simplification of the inverse transformations is 
possible in orbital analysis by associating the elliptic coordinates 
with the orbital elements c 

Consider the equation of an ellipse in the form 

* 2 + ^ = i 

dr cosh u d sinh u 


(9) 


Equation (9) is often used in the derivation of Equations (8 ) ^ The 
standard equation of an ellipse is given by 


% + - 1 CIO) 

a t> 


where a is the semimajor axis and b is the semiminor axis. By 
comparing Equations (9) and (10) , it is evident that 
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a = d cosh u 
b = d- sinh u 


Cll) 

( 12 ) 


Consider now the expression for the eccentricity in the form 




2 , 2 - 
a - b 


(13) 


Substituting Equations (11) and (12) into Equation (13) yields 


e = 


d“ cosh 2 u - d 2 sinh 2 u* 


2 — 

d cosh z u 


(14) 


which reduces to 


_ j- __ ! ^ 

■cosh 2 u-l 


■ — i 


or 


cosh u = 


(16) 


Substituting Equation (16) into -Equation (11) yields 


d = ae 


(17) 


Note that if u is chosen as one of the coordinates, 
then from Equation (16) u is given by 


u = cosh 


( 18 ) 


which is singular for circular orbits (zero eccentricity) . However, 



by defining the coordinate as d cosh u, this singularity is elim- 
inated since d cosh u = a as shorn in Equation (11) . 

Now, recall that the parametric equations. of an elliptic 
orbit are given by 

x = a cos E y = b sin E 

Solving Equations (11) and (12) for d and then substituting these 
equations into the respective terms of Equation (8) yields 

.X = a cos v = b sin v 

e e 

A comparison of Equations (19) and (20) demands that 


Hius, since d cosh u = a and v = E, coordinate sets 1 and 2 may 
be given by 


“1 = a 



a 3 


= z 


and 


B^ = a $2 ~ cos E = z 

Note that the coordinates as given in Equations (22) and (23) re- 
quire that the xy and uv planes coincide with the orbit plane. 
Thus, since and are initially zero and since the motion 



2 i 


is constrained to the orbit plane, “ 3 >“ 3 >£ 3 >63 mu5t ^ en a ^ 
necessarily be zero. In addition, since a-^ and pj are deter- 
mined from the semimajor axis which is constant for unperturbed 
orbits, the terms are a l so zero> As a result of 

these simplifications, the only coordinates with nonzero accele- 
rations are and If the time derivatives of these two 

coordinates are expressed in terms of the eccentric anomaly, the 
accelerations for these coordinates may be given by 


a-, = 


-ye sin E 


2 'SfT. c', 3 

a (1 - e cos E) 


(24) 


and 


V* i->j 


a^(l - e cos E)~* 


(25) 


Thus, the number of coordinates to be integrated has been reduced 
from three to one. 

Since the elliptic coordinates are referenced to the cen- 
ter of the ellipse (Figure 5) and since elliptic orbit motion is 
referenced to the focus of the ellipse (Figure 1) , the origin of 
the elliptic coordinate system must be transferred to the focus. 

This is accomplished by translating the x-axis according to 

x = a-^ cos a 2 - ae (26) 


or 
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x = 3 ^ " ae ( 27 ) 

Note that since coordinate set 3 is referenced to the equa- 
torial plane, the simplifications as given in Equations (11) and (12) 
are not applicapable to these coordinates. This requires the coor- 
dinates to be determined by solving Equation (8) , which is quite 
cumbersome. Furthermore, when referenced to the equatorial plane, 
all three coordinates must be integrated instead of one as in sets 
1 and 2. For these reasons coordinate set 3 is discarded in this 
study. 


11,3,6 Ellipsoidal and Conical Coordinates 
The ellipsoidal and conical coordinates were discarded 
necause ox pj ctt ' amhi gurries arm mg j a me coord in^ ,p tran ^Tryrmji - ' 
tions. The transformations from the rectangular Cartesian to the 
ellipsoidal and conical coordinates were also quite cumbersome. 

11. 4 Equations of Motion 

The two -body equations of motion for the coordinate systems 
investigated in this study are presented in Appendix B. The equa- 
tions of motion for the rectangular and the spherical coordinates 

were obtained from Reference 1. The equations of motion for circular 

\ 

and elliptic coordinates were obtained using Lagrange's equation. 

11. 5 Computational Procedure 

The computer programs used in this study were written in 
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Fortran IV for the CDC 6600 computer (listings of the program are 
available upon request) . Since this computer has a 60-bit word 
length and since single and double precision Univac 1108 results 
were available as checks, it was felt that double -precision arith- 
metic would be unnecessary. In order to verify this assumption 
some of the work was carried out in both single and double pre- 
cision. No appreciable difference in accuracy was found. How- 
ever, the double precision results did require a longer computa- 
tional time as would be expected. The computational algorithm as 
* 

rused in this investigation is as follows: 

1. Prior to any major calculation the transformation 
from the orbital plane to the reference coordinate system is 

'"gsneT.'red. 

2. The initial reference state in the orbit plane and the 
transformation to the reference coordinate system are then generated. 

3. The initial conditions for each coordinate system are 

i 

then set up using the necessary states formed in step two. 

4. The equations of motion for each set of coordinates for 
the specified integration period are integrated numerically. At 

. specified intervals the appropriate error checks are performed. The 
numerical integration routine (which was obtained from Fowler and 
Lastman - Reference 3) uses a fourth order Runge-Kutta as a starter 
and then shifts to a fifth order Adams predictor-corrector. As used 
in this investigation, the entire integration procedure was carried 
out in single precision. 



CHAPTER III 


RESULTS AND DISCUSSION 


The results of this investigation are presented in the forms 
of graphs (Figures 6 through 19) and tables (Tables 6 through 11) 
showing the error in position and velocity as a function of orbit 
time for various combinations of step size and orbit shape. Since 
the errors in state range over several orders of magnitude, each 
coordinate error curve is scaled (if necessary) so that all the 
coordinate error curves for one step size and orbit shape combina- 
tion can be presented on one graph. If scaling is necessary, the 
respective scale factor is indicated on each graph in paran thesis. 

iJuxoao o-jljl lautulo alo OjuiJo jlij. uiuui Cv 

obtain the correct error values for the scaled system, simply 
divide the ordinate value by the indicated scale factor. An 
example of the raw computer plot output is shown in Figures 6 and 
13 with Figure 6 showing the position error and Figure 13 the veloc- 
ity error. The large oscillations of the coordinate error curves 
cause the curves to overlap. Therefore they do not emphasize the 
item of primary concern; i.e., the size of the maximum error in the 
vicinity of a giiren point in time. This is most easily shown by 
plotting maximum error envelopes for the various error curves. The 
envelopes of the curves shown in Figures 6 and 13 are shorn in Fig- 
ures 7 and 14. Since the graphs showing the maximum error curve 
envelopes are interpreted more easily than plots showing the raw 
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ELLIPTIC COORDINATES, SET 1, TERMINAL 
TIME VALUES FOR CIRCULAR ORBITS 

Step Size 

Position Error 

Velocity Error 

(seconds). 

(feet) 

(feet/second) 

15 

0.07284 

0.00008698 

30 

0.02406 

0.00002873 

60 

0.02120 

0.00002531 

120 

0.01546 

0.u00u1846 

150 

0.00309 

0.00000369 

240 

0.00678 

0.00000809 


TABLE 6 




ELLIPTIC COORDINATES TERMINAL TIME COMPARISONS 


Position Error 
(feet) 


Set 1 


11.3869 


Set 2 

Set 1 

11.2482 

0.0004 

12.0965 

0.0123 

175.8160 

0.0032 

172.6260 

0.1810 

i 

2oyl.4/uu 

u.u4/i 

1986.0400 

0.3484 

38759.4000 

0.6501 

9142.8400 

10.5715 


Velocity Error 
(feet/second) 


Set 2 


TABLE 7 
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ERROR COMPARISON 15 SEC STEP SIZE 

A/P 

R 

C 

S 

E2 

El 

Terminal Time 

Position Error (feet) 

1 

100.910 

6.5387 

2S.3569 

12.4850 

0.0728 

2 

77.029 

15.6925 

33.0283 

11.2482 

0.3196 

5 

. 15.298 

28.3731 

36.9600 

9.9504 

1.0695 

10 

83.986 

44.7332 

42.9504 

8.7478 

2.9602 

20 

271.134 

97.8794 

52.8845 

12.0965 

11.3869 



Maximum Position Error (feet) 


1 

100.910 

18.0245 

30.8165 

17.9850 

0.0728 

2 

77.248 

19.7348 

33.2170 

13.8808 

0.3196 

5 

36.105 

33.4939 

38.8750 

14.6235 

1.0695 

10 

125.670 

58.1618 

44.1784 

90.0689 

3.2218 

20 

271.134 

97.8794 

52.8845 

86.6879 

11.3869 

\ 

T n Vol nri tv Prrnr ( -Ft* / cor 

*/ " . * 

" j 

1 

0.1205 

0.0212 

0.0367 

0.0149 

0.00009 

2 

0.0850 

0.0150 

0.0367 

0.0131 

0.00037 

5 

0.0145 

0.0322 

0.0431 

0.0112 

0.00121 

10 

0.0846 

0.0465 

0.0455 

0.0090 

0.00305 

20 

0.2712 

0.1002 

0.0591 

0.0132 

0.01230 

Maximum Velocity Error (ft/sec) 

1 

0.1205 

0.0214 

0.0367 

1.3083 

0.00009 

2 

0.0850 

0.0222 

0.0379 

2.6678 

0.00038 

5 

0.0331 

0.0331 

0.0431 

2.4131 

0.00124 

10 

0.1218 

0.0630 

0.0519 

51.5267 

0.00366 

20 

0.2712 

0.1002 

0.0591 

14.6945 

0.01230 


TABLE 8 
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ERROR COMPARISON FOR 30 SEC STEP SIZE 


A/P 

R 

C S E2 El 

Terminal Time Position Error (feet) 

1 

3415.7800 

417.1100 

494.0320 

172.6260 

0.02406 

2 

2839.1400 

385.8500 

565.8650 

175.8160 

2.73115 

5 

1201.2200 

118.0900 

626.0040 

159.0570 

14.06290 

10 

372.5100 

400.5900 

709.5340 

134.9860 

42.59470 

20 

2920.4300 

1205.3400 

846.8990 

172.6260 

167.52200 


Maximum Position Error (feet) 

1 

3415.7800 

487.1000 

529.9230 

389.4900 

0.02400 

2 

2839.1400 

407.0700 

566.7660 

326.4700 

2.73900 

5 

1651.3300 

246.5800 

652.4170 

211.4900 ! 

14.06300 

10 

1060.4600 

589.0800 

733.2990 

1039.2900 ; 

46.01300 

20 

2920.4300 

1203.3400 

846.8990 

1373.5400 

167.52200 


7'erminal Time Velocity Error (ft/sec) 

! 

1 

4.0/85 ; 

0.5864 

i 

0.6320 

0.2387 

0.00003 

2 

3.2226 

0.3842 

0.6298 

0.2054 

0.00317 

5 

1.4538 

0.2175 

0.7287 

0.1807 

0.01595 

10 

0.3676 

0.4303 

0.7512 

0.1387 

0.04382 

20 

2.8086 

, 1.2096 

0.9449 

0.1920 

0.18098 


Maximum Velocity Error (£t/sec) 

1 

4.0785 

0.5864 

0.6320 

42.7050 

0.00003 

2 

3.222$ 

0.4414 

0.6473 

119.1500 

0.00325 

5 

1.5252 

0.2323 

0.7287 

30.6480 - 

0.01635 

10 

0.9935 

0.5687 

0.8596 

581.5750 

0.05274 

20 

2.8170 

1.2096 

0.9449 

230.1740 

0.18098 


TABLE 9 
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ERROR COMPARISON FOR 60 SEC STEP SIZE 



C S E2 

1 

Terminal Time Position-Terror (feet) 


108908.000 

93982.800 

60304.500 

33929.000 

2060.600 


25043.7000 

10096.2000 

3222.5300 

21903.1000 

11009.7000 

2691.4200 

12927.8000 

11728.4000 

2560.5300 

5448.4000 

12519.4000 

1998.5500 

6238.0000 

13756.1000 

1986.0400 


0.02000 

40.57000 

200.91000 

582.38000 

2173.85000 



Maximum Position Error (feet) 


1 

108908.000 : 

25131.1000 

10589.4000 

9899.8000 

0.02000 

2 

93982.800 

21903.1000 

11009.7000 

13213.5000 

40.57000 

5 

641^-3.200 

14619.0000 

12014.7000 

54164.9000 

200.91000 

10 

40482.500 

8070.5000 

13095.0000 

20967.4000 

635.33000 

20 

22929.800 

6493.0000 

13756.1000 

20334.2000 

2173.85000 

Terminal Time Velocity Error (ft/sec) 

1 

130.034 

30.2844 

12.5889 

3.8593 

0.00003 

2 

107.814 

24.7957 

12.3223 

3.2326 

0.04707 

5 . 

69.674 

15.4015 

13.5759 

2.9181 

0.22788 

10 

35.999 

6.4585 

13.2545 

2.0487 

0.59914 

20 

5.972 

5.1896 

15.3152 . 

2.5780 ! 

2.34835 

Maximum Velocity Error (ft/sec) 

1 

130.034 

. 

30.28440 

12.5889 

1590.2700 

0.00003 

2 

107.814 

24.79570 

12.5631 

7419.3900 

0.04827 

5 

69.674 

15.40150 

13,5759 

216.8600 

0.23411 

10 

36.344 

6.45846 

15.2794 

5793.7300 

0.72210 

20 

22.150 

5.79800 

15.2794 

2920.7800 

2.34835 


TABLE 10 























46 



ERROR COMPARISON 

FOR J 20 SEC STEP SIZE 



R 

C 

S 

E2 

El 

A/P 








Terminal Time Position Error (feet) 


1 

3064172.00 

1278348.000 

379119.000 

53863.1000 

0.02000 

2 

2708974.00 

1127995.000 

364646.000 

38759.4000 

560.32600 

5 

2059712.00 

856197.000 

339097.000 

42396.6000 

2205.58000 

10 

1586616.00 

628725.000 

301231.000 

24741.3000 

4765.60000 

20 

1592896.00 

547540.000 

231853.000 

9142.8400 

9792.56000 

Maximum Position Error (feet) 

1 

3064172.00 

1278348.000 

380589.000 

594284.0000 

0.02000 

2 

2708974.00 

1127995.000 

364646.000 

243501.0000 

564.63000 

s • 

2059712.00 

856197.000 

339097.000 

197259.0000 

2205.58000 

10 

1619416.00 

649060.000 

322389.000 

175354.0000 

5268.40000 

20 

1592896.00 

547540.000 

231853.000 

195278.3000 

12389.52000 



Terminal Time Velocity Error (ft/sec) 


i 

ODDO.iU 

jSz/.u4u 

.. .. 
4^0. 1U4 

o5.u/zi 

. ! 

U.UUUU.£ 

2 

3116.62 

1298.000 

416.088 

55.5405 

0.65011 

5 

2344.99 

975.208 

386.801 

49.0478 

2.50149 

10 

1670.24 

661.111 

317.481 

25.5778 

4.90396 

20 

1809.25 

616.306 

256.427 

59.9667 

10.57150 

Maximum Velocity Error (ft/soc) 

1 

3656.20 

1527.040 

456.104 

12324.3000 

0.00002 

2 

3116.62 

1298.670 

418.667 

9012.9000 

0.66885 

5 

2344.99 

975.210 

386.801 

15553.2000 

2.59316 

10 

1896.08 

761.480 

372.850 

41486.6000 

5.98710 

20 

1809.25 

616.310 

256.427 

12112.5000 

10.57150 


TABLE 11 
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data, only the plots showing the maximum error envelopes will be 
shown for the remaining data. However in discussing the results, 
reference will still be made to the oscillation in the error curves. 
The amplitudes of these periodic deviations increase for both step 
size increases (constant eccentricity) and eccentricity increases 
(constant step size). 

III.l Effect of Step Size on Error Propagation 

In analyzing the data generated in this investigation for 
step size effects (Figures 7, 8, and 9), one characteristic imme- 
diately noticeable is that the error plots for each coordinate 
system are oscillatory in nature with a period corresponding to 
that of the orbit as is shown in Figure 6. Closer examination 
reveals that the amplitudes of oscillation are dependent on the 
integration step size for constant orbit shape. This dependence 
is such that as the step size increases, the oscillation amplitude 
also increases. With but one exception, all of the coordinate 
error curves exhibit this type of periodic deviations superimposed 
upon nearly linear increasing functions of orbit time. The rectangu- 
lar coordinate errors increase in a nonlinear fashion with large 
periodic deviations. In the following paragraphs the relative 
influence of step size on each coordinate system is discussed. 

The nonlinear characteristic of the average rectangular 
coordinate error curves indicate that this system is indeed sensi- 
tive to integration step size. The rectangular coordinate error 
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usually increases slowly with orbit time for the smaller step sizes 
(constant eccentricity) as shown in Figure 6 or 7. As the integra- 
tion interval increases, the rectangular coordinate error reaches 
a local maximum (Figure 8) and then declines to a local minimum. 

Upon passing through this local minimum, the state deviations in- 
crease almost exponentially (Figure 9) . Increasing the step size 
further causes the local maximum and minimum to be reached much 
earlier in time (Figure 10) thus allowing more time for terminal 
error buildup. This trend is not as apparent in the other integra- 
tion error curves, although the errors in the circular coordinate 
system exhibit a greater step size dependence than do the spherical 
coordinate errors. 

In comparing the circular coordinate errors to the spherical 
coordinate errors , it is apparent that the errors which occur in the 
spherical system are generally larger than those which occur in the 
circular system for the IS and 30 second step sizes and small eccen- 
tricities as shown in Figures 8, 11, and 12. However, for the larger 
eccentricities the spherical error curve approximates the average 
circular error curve (Figure 6). For these cases, the circular error 
curve exhibits large deviations about the average curve , For the 
larger step sizes, 60 and 120 seconds, the spherical deviations are 
larger than the circular errors early in the accumulated integration 
time. Hi is occurs prior to the buildup of truncation error result- 
ing from too large a step size (Figure 10) . As the effects of trun- 
cation error increase, the circular errors become larger than the 
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spherical coordinate errors as is expected. This effect leads to 
the conclusion that circular coordinate errors .are more sensitive 
to step size than either the spherical or elliptic coordinate 
errors . 

The coordinate errors in the elliptic coordinates, sets 
1 and 2 (as defined in Equations 5 and 6 , respectively) , generally 
increase only slightly as the step size increases are made. Coor- 
dinate set 2 exhibits a greater error sensitivity (error increase) 
to step size for all orbit shapes than does coordinate set 1. 

In addition, the 'scattering of a few data points (the few extremely 
large deviations from the average error curve - see Figures 6 and 13) 
in set 2 increases in magnitude as the step size increases. The few 
scattered points were not considered in tne determination of the 
maximum error envelopes for set 2. Coordinate set 1 appears to be 
only slightly sensitive to step size with the errors increasing more 
for the larger eccentricities. Furthermore, coordinate set 1 does 
not exhibit any data scatter. An exception to the error increase 
with step size increase is noted for circular orbits. As the step 
size increases, the state deviations decrease. Further investigation 
shows that the errors in the state decrease to a minimum between the 
step size of 120 seconds and 240 seconds. Note that the terminal 
time comparisons given in Table 6 are valid, since for circular 
orbits, the errors of the elliptic cylinder coordinates, set 1, in- 
crease without oscillation. This indicates that the "best" (most 
accurate) step size for circular orbits exists in the step size range 



of 120 to 240 seconds for the CDC 6600 Computer. The decrease in 
state errors for step size increases is evidently due to a reduction 
in round-off error resulting from fewer integrations. The state error 
increases for the larger step sizes must be due then to truncation 
error. Note that none of the other coordinate systems investigated 
indicates the existence of a "best" step size in the range con- 
sidered (15 seconds to 120 seconds) regardless of orbit shape. It 
follows then that the "best" step size for any of the other coordinate 
systems considered (if it exists) must be less than 15 seconds. This 
indicates that the "best" step size (maximum -accuracy) depends upon 
the word length of the computer involved with the larger word lengths 
allowing the most accuracy in the integrated states. This hypothesis 
is further supnorted bv that fgrt that Cerhep an a Tewgilen fe ,j nd 
that for circular orbits a best step size of 30 seconds existed for 
the spherical coordinate system for the Univac 1108. 

In general then, as the integration step size increases, 
both the position and velocity errors increase for all combinations 
of coordinate system, orbit shape and step size. Hie lone exception 
is for circular orbits using set 1 of the elliptic coordinates as 
previously mentioned. 

III. 2 Effect of Orbit Shape on Error Propogation 

In examining the plotted data for eccentricity effects, it 
is observed than the error curves, in addition to being step size 
dependent, are also eccentricity dependent. Hie state deviations 
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display small amplitudes for small eccentricities (constant step size) 

and approximate a smooth curve for circular orbits. As the eccentricity 
« ‘" Ir 

increases, the amplitudes also increase with the rectangular coordinate 
system exhibiting the largest increase. 

The amplitudes of the oscillating rectangular coordinate 
errors are highly eccentricity dependent increasing significantly as 
the eccentricity increases. In addition to influencing the error 
curve amplitudes, the eccentricity increases affect the average error 
curve in a reverse manner to that of the step size variations. In 
this case, as the eccentricity increases (constant step size) , the 
local maximum and minimum begin to appear later in orbit time (Fig- 
ures 8, 11, and 12), However, an increase in the amplitude of the 

^ A ^.-^4 ^4- 4 ~ ^ 4-t 4.^. ^ * i rnl. • _ 4 . 1 ^ 
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ably results from the fact that the coordinates vary more slowly for 
the larger orbits than for the smaller orbits. 

A close examination of the effect of orbit shape upon the 
spherical and circular coordinate error curves indicates that both 
of these systems are affected in a similar manner. This conclusion 
is based on the fact that early in the orbit time, when step size 
influence is minimal, the spherical coordinate errors are usually 
larger than the circular coordinate errors. As truncation error 
influence (due to large step size) increases, the circular errors 
become larger than the corresponding spherical errors (Figure 10) . 

This trend continues for nearly all eccentricities. Thus, it is 
concluded that eccentricity does not appear to influence the state 
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deviations in the spherical and circular coordinate systems to a 
high degree. This same conclusion can be extended to the elliptic 
coordinates . 

The elliptic coordinates, sets 1 and 2, both appear to be 
the coordinate systems the least influenced by eccentricity changes. 
Set 2 continues to exhibit scattered data points in the vicinity of 
apogee. Coordinate set 1 appears to be slightly more sensitive to 
eccentricity due to the fact that the magnitude of the terminal 
position errors increases significantly as eccentricity increases. 
Table 7 compares state errors at the terminal time for both systems 
indicating to a degree the sensitivity of each to eccentricity. 
However, note that the terminal time comparisons are not completely 
valid d 11 ? to the n sci il at?nc nature nf the error curves. Note also 
that although both coordinate sets have approximately the same 
error, set 1 would still be more desirable computationally because 
its error envelope is much smaller than that for set 2. 

It may thus he concluded that all the coordinate error 
curves are influenced by eccentricity. This influence though , is 
not as strong as that due to step size. 

I II, 3 Effect of Coordinate System on Error Propagation 

The effect of coordinate system alone upon the errors in 
state is indeed difficult to isolate due to .the interaction of the 
orbit shape and the integration step size. In the paragraphs to 
follow, each coordinate system will be separately discussed, and the 
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respective advantages or disadvantages given. 

line rectangular coordinate system is highly sensitive to 
step size and eccentricity. Generally speaking "the error magni- 
tudes become quite large and dominate the error magnitudes in all 
the other coordinate systems considered. This is evident in Tables 
8, 9, 10, and 11 which compare terminal time error with the largest 
error occuring during the integration interval. Due to the large 
errors the rectangular coordinate system is not recommended for 
numerical studies. 

Hie circular coordinate system is also sensitive to in- 
creases in step size and orbit shape as indicated in Tables 8, 9, 

10, and 11. Although the circular coordinate errors are small in 
some instances, the large periodic oscillating error curves generally 
mean large state errors. A comparison of the circular coordinate errors 
with those obtained using spherical coordinates indicates that both 
systems have potential uses depending upon the type of orbit to be 
investigated and the step size to be taken. 

The spherical coordinate state deviations are quite stable 
and do not vary drastically for step size changes (see Tables 8, 9, 

10, and 11). Orbit shape appears to have little effect on coordinate 
errors in this system. The stability of this system indicates that 
it may be very useful for investigating perturbed orbits. 

The elliptic coordinate systems produce the most accurate 
states (Tables 8, 9, 10, and 11). The scattering of isolated data 
points in coordinate set 2 is highly undesirable. Coordinate set 1, 
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.however, is generally quite accurate . -In addition this system also 
appears to be a desirable one in which to investigate perturbed 
motion. One possible objection to this system Us used here, is that 
the system is rectified to the orbit plane. Tnis requires a coordi- 
nate transformation to obtain the inertial states. However, the 
accuracy gained more than outweighs this disadvantage. It thus 
appears that the elliptic cylinder coordinates, set 1, offer superior 
computational accuracy and exhibit computational stability over long 
time intervals even for large step sizes. 



CHAPTER IV 


CONCLUSIONS AND RECOMMENDATIONS 


IV. 1 Summary ^ 

The sensitivity of two body states obtained through numerical 
integration to coordinate system and integration step size is investi- 
gated for various orbit types. This is accomplished by numerically 
integrating the equations of two-body motion for the various coordinate 
systems and then comparing resulting states to the states generated 
through the two-body solution. By holding one parameter constant and 
varying the other two, the effect of each on the integrated states can 
be determined. It is observed that for all coordinate systems con- 
sidered in this investigation, step size has a greater influence on 
the error than- the orbit shape. This may be seen by comparing Figures 
7, 3, and 9 to Figures 8, 11, and 12. Also, it can be observed that 
the position and velocity error curves behave in a similar manner. 

, This is best seen by comparing Figure to Figure . Finally, any 
operation designed to reduce the errors in position also reduces the 
errors in velocity as is evidenced by the corresponding decreases in 
position and velocity errors of the elliptic cylinder coordinates 
(Tables 8 , 9 , 10 , and 11) . 

IV . 2 Conclus ions 

Based on the data generated in this investigation the follow- 
ing conclusions have been reached. Some of these conclusions are made 
with reservations since additional investigation is necessary before 
final conclusions are possible. 
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1. Elliptic coordinates, set 1, exhibit superior computa- 
tional accuracy when compared to the other coordinate systems. This 
system remains relatively insensitive to orbit shape and integration 
step size as far as error increases are concerned. The error curve 
increases in a linear manner, and the periodic variations superimposed 
on the ave rage error curve are generally quite small. The amplitude 
of the periodic deviations increases only slightly, thereby providing 
a small error envelope. 

2. Elliptic coordinates, set 1, allow much larger step 
sizes to be taken before similar error magnitudes occur. These large 
step sizes represent a cost decrease by reducing the running time. 

3. Elliptic coordinates, set 2, also produce error curves 

t *1% /v« h i n r»Y« 1 1 T ir\^Tf\-v+ i caI o + q/ 1 qcotfAror 
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error points cause the error envelope to become quite large. This 
is the only unattractive feature of this system. Note, however, 
that even without the isolated scatter of error points, set 1 is 
still more accurate than set 2. 

4. Given the option of choosing from a set of coordinate 
systems, one containing the angle itself and the other containing 
the cosine of the angle as one of the coordinates, the set containing 
the angle should be chosen. Integration of the angle is inherently 
more accurate than integration of the cosine since the angle increases 
monotonically while the cosine is oscillatory. This is verified by 
comparing the results of elliptic coordinates, set 1, to elliptic 
coordinates, set 2. Both sets are identical in all respects except 
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that in set 1 the angle is integrated while in set 2 the cosine of. 
the angle is integrated. Thus, any difference in error magnitudes 
between the two systems must result from the integration of the 
cosine . 

5. The spherical coordinate error curves are well behaved 
and increase linearly with time. Periodic deviations of almost con- 
stant amplitude are superimposed on the average error curve for all 
combinations of orbit shape and step size. Hie small periodic devia- 
tions about the linear average error curve is a desirable feature of 
this system because the small deviations keep the error envelope 
small thereby allowing the actual error curve to be approximated by 
the average error curve. 

Til a r/'lll OV j-rt r 4-rt 
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- ing to step size and orbit shape. Depending upon the step size and 
orbit shape, the coordinate errors either remain small or increase 
rapidly with large periodic deviations which cause the error envelope 
to become quite large. This type of error propagation is the un- 
desirable feature associated with using this coordinate system. 

7. For small step sizes (15 and 30 seconds) and small orbits 
(A/P = 1, 2, 5), the circular coordinate errors are usually smaller 
than the spherical state deviations. For large orbits (A/P = 5, 10, 

20) and both large and small step sizes, the spherical coordinate 
coordinate errors are generally smaller than the circular state errors. 

8. Rectangular coordinate error curves are greatly influenced 
'by step size and orbit shape and usually produce' extremely large errors. 
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These large errors are the results of either the large periodic varia- 
tions (large error envelope) superimposed on the average error curve 
or the nonlinear (exponential type) increase of* the average error 
curve . 

9. The "best" step sizes (those step sizes for which the 
integrated states are the most accurate) appear to be dependent upon 
coordinate system and computer word length. 

IV. 3 Recommendations 

The following recommendations for further study are made: 

1. Similar data should be generated on other computers 
in order to determine how step size as a function of word length 
effects the accuracy of the integrated states. 

2. A thorough study in which other coordinate systems 
than those considered in this -study are included should be made for 
a wider range of integration step sizes and orbit shapes. The com- 
putation involved should be done on several different computers 
using both single and double -precision arithmetic in order to deter- 
mine the effect machine and arithmetic type on the accuracy of the 
integrated states. 
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APPENDIX A 


The transformations between coordinate systems as used in 
this study are presented in the following sections. 


A.l Transformations Between Spherical and Rectangular Coordinates 
Spherical to Rectangular 

X = p cos a cos t Y =p cos a sin t Z = p sin a (1) 

and by differentiation 

X = p cos a cos t - p d sin a cos x - p x cos a sin x 
Y = p cos a sin t - p d sin a sin x + p x cos a cos x (2) 

Z = p sin o + p b cos a 

Rectangular to Spherical 
Again, fr~ gconexry 
p = (X 2 + Y 2 + Z 2 ) h 

x = tan _1 (Y/X) (3) 

a = tan~ 1 (Z/(X 2 + Y 2 ) h ) 

and by differentiation 

. . XX + YY * ZZ 

P + 7 7 7 U 

Or + y z + i 1 ) 2 


. _ XY - YX 
x 7 7 

(X^ + Y j 

• = CX 2 + Y 2 )Z - Z(XX + YY) 
(X 2 + Y Z )' 2 (X 2 + Y 2 + Z Z ) 


( 4 ) 
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A. 2 Transformations Between Circular and Rectangular Coordinates 
Circular to Rectangular 
From Figure 4,, 

X = r cos 3 Y = r sin 02=2 (5) 


and by differentiation, the velocity transformations become 

X - r cos 0 - r 3 sin e Y = r sin e + r e cos e 
Z = 2 


( 6 ) 


Rectangular to Circular 
Again, from geometry, 

r = (X 2 + 8 = tan -1 (Y/X) z = Z 

and by dixxereiu.iat.xoju the velocity TiXmaxuliiiatxOUa bcCwiiiO 


XX + YY 

r = „ 
(T + Y 



XY - YX 

~2 ~T 
X + Y 



C7) 


( 8 ) 


A. 3 Transformations Between Elliptic Coordinates, Set 1, and 
Rectangular Coordinates 

Elliptic Coordinates, Set 1, To Rectangular 

2 % 

X = cos ^2 - ae Y = oj, (1 - e ) 2 cos 

(9) 

Z = o 3 


and by differentiation 
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V 

a = cos c *2 - a i a 2 sin a 2 

, . 7 U • 

Y = a( 1 - e' ) sin “2 + a i a 2^ " e ^ cos a 2 (10) 



Rectangular to Elliptic Coordinates, Set 1 

= a 0.2 - E .a^ = 0 (11) 

and by differentiation 

a-^ = 0 “ E = 0 (12) 


A. 4 Transformations Between Elliptic Coordinates, Set 2 , 
and Rectangular Coordinates 


bilipric uoor ujuiates } Set 2 > to i\ectaiigu.i-ia.i‘ 


X = 


e i 6 2 


- a e 


Y = B, (1 - 


7 l- 

0*(1 




Z = 3, 


(13) 


and by differentiation 


X ~ 3^3 2 ^1^2 


z = 8, 


Y = 


^[(1 - e 2 )(l - 3 2 )]^ Y - S^3 2 6 2 (l - e 2 ) 


( 14 ) 


Note that Y does not exist for Y equal zero and must be determined 
by another method (h = r x v) . Note also that the proper sign for Y 
is ambiguous and must be determined from considerations of X and X. 
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Rectangular to Elliptic Coordinates , Set 2 

&1 = z B 2 - cos E = 0 (15) 

and by differentiation 

3-^=0 B 2 = "E sin E = 0 (16) 


A. 5 Transformation From the Orbital Plane to the Reference 
Coordinate System 

The transformation from orbital elements to the inertial 
rectangular coordinate system is well known and is derived by various 
authors. The transformation will not be derived here but is given as 

j c- — 1 

UL/L.UXUWVJ i 7 wilt iVt# t ^ 


” x“ 


a ll a 12 a 13 


~ x" 


Y 


a 21 a 22 a 23 


y 

(17) 

Z 

Inertial 

a 31 a 32 a 33_ 


z 

Orbital 





Plane 


where 

= cos co cos ft - cos i sin ft sin co 

a 12 ~ ' s ^ 11 40 co5 n ” cos * s * n ® cos 03 
a 13 ~ s * a a s ^ n 

a 2 ^ = cos to sin ft + cos i cos ft sin co 
a 22 = 


sin to sin ft + cos i cos ft cos w 
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a 23 = 

- sxn 

1 cos 

a 31 = 

sin 

U) 

sin i 

a 32 = 

cos 

Cl) 

sin i 

a 33 = 

cos 

i 





the angles i, n and u are the orbit inclination, the longitude of 
the ascending note and the argument of perigee respectively, The 
coordinates x, y and z are given by 

x = a cos E - a e y - b sin E z = 0 (18) 
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APPENDIX B 
EQUATIONS OF MOTION 

The equations of motion as used in this study are presented 
in the following sections. 


B.l Rectangular Coordinates 



where r - (X^ + 


Cl) 


B . 2 Spherical Coordinates 

n = — ii- + cos^ a + ncf^ 

V 

t = 2ra tan o - 2pi/p (2) 


-2 

a .= -t cos o sin a - 2pa/p 


B.3 Circular Coordinates 


•• pu *2 

P ' ' ITT 7y7Z P 6 


2p0 v - "P 2 


( 3 ) 
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B.4 Elliptic Coordinates, 

a l = 0 a 2 

B.5 Elliptic Coordinates > 
3-^ = 0 g 2 


Set 1 


“pe sin v 

3,, 3 

a (1 - e cos v) 



( 4 ) 


Set 2 


u(e - cos v) 
a (1 - e cos v) 



C5) 
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